﻿ Prelucrarea semnalelor Facultatea de Electronica si Telecomunicatii, UPT http://shannon etc upt ro/teaching/ps/ 1 Transformarea Laplace http://shannon etc upt ro/teaching/ps/Cap7 pdf Pierre Simon Laplace 2 Regim permanent – Transformata Fourier Regim tranzitoriu –Transformata Laplace 1 •Transformarea Laplace bilaterala j st 11 DC,tsXdsesXtxL  2j j  st ;jdtetxsstxXL ;s ,    •Transformarea Laplace unilaterala   j 1 s t  s t                Re sx t X s e ds s x t e dt0;   02jj   •Relatia dintre cele doua: stxstxstxLLL uu3 Trasformata Laplace generalizează transformata Fourier în tot planul complex Pentru variabila complexă se utilizează notaţia: (1 1)    ,     js  j =cst  4 2 Pentru =0 transformata Laplace este spectrul semnalului adică transformata Fourier j =cst  5 1 Răspunsul unui sistem liniar şi invariant în timp la o exponenţială cu exponent complex Semnalul exponentială complexă, cu exponent complex:  tjtts (1 2)     ,000eeetx se aplică la intrarea unui SLIT caracterizat prin răspunsul la impuls h(t) La ieşirea sistemului se obţine semnalul:      00jt y t h t x t h x t d h e d                          0 0 0 00j t jsts e h e d e h e d               0ss In ipoteza că este convergentă, se notează cu H(s) integrala:  js deehdehsH  (1 3)  - depinde numai de h(), 6 - în anumite condiţii caracterizează complet SLIT 3 Răspunsul y(t) al SLIT : 0st (1 4)     00y t e H s x t H s st care ne arată că : - 0eeste o funcţie proprie a SLIT, - H(s) este valoarea proprie asociată ei 0 Daca semnalul de intrare este o combinaţie liniară de forma: ts k ,ectx(1 5) k k rezultă : ts k(1 6)      , kkesHcty k operatorul S care modelează sistemul fiind liniar Se poate vedea că o funcţie H(s) definită prin expresia (1 3) permite calculul răspunsului unui SLIT la care semnalul de intrare poate fi pus sub forma (1 5) 7 2 Transformarea Laplace bilaterală Prin definiţie transformata Laplace bilaterală a unui semnal x(t) este:  st(1 7) LX  x t s s x t e dt s j; ; ,                     Legătura între transformarea Laplace bilaterala şi transformarea Fourier:  t j t t    (1 8) x t j x t e e dt x t e      LF       Transformata Laplace bilaterală a unui semnal x(t) este transformata t Fourier a semnalului x(t)e- Daca =0, se obţine:                     XtxjjtxFXL(1 9)    Transformata Laplace bilaterală evaluata pe axa imaginara jω (σ=0) este 8 egala cu transformata Fourier a aceluiasi semnal 4 Notaţii: stxXL În loc de a scrie , se scrie (din raţiuni tipografice) prin urmare:  sXtxL şi txjXjtxFL în continuare, vom nota cu X(s) transformata Laplace bilaterală, tolerând două notaţii pentru transformata Fourier şi anume X() când este definită direct şi X(j) când rezultă prin particularizarea transformatei Laplace bilaterale:        XjXsX(1 10)     js  9 Exemple:  t     0tetx 1 ) Considerăm semnalul cauzal  Transformata 1  txF       sa Fourier există şi este:    j 0 Transformata Laplace bilaterală a semnalului x(t) este:           00ttst j t L    x t s e t e dt e e dt          0 Integrala este convergentă numai pentru +>0, deci  > -, 00 care se mai scrie şi Re{s}>-  Avem deci: 0 11 (1 11)        x t s s0; Re L       00js Se poate constata că, dacă - -2, atunci fiecare termen al relaţiei (1 16) va fi inversat într-un semnal cauzal: 21t 2 : et     ;  Fig 1 5 Cele 3 DC posibile s2 pentru X(s)=1/[(s+2)(s+3)] 31t -3 : et     s3 Rezultă: 2311tt 2 : e e t          (1 17) 20 s 2 s 3 10 ii) Dacă DC este 2,  -2, şi zerouri a unei situat la dreapta verticalei ce trece prin poli transformate Laplace Axa j este inclusă în DC există şi transformata Fourier a semnalului:   00sjss F   xt                    2 1 21p p p pssssssjj sj Transformarea Fourier este particularizarea lui X(s) for s=jω Ea poate fi Pentru fiecare valoare, punctul A ia o calculata pentru diferite valori ale ω24 pozitie diferita pe axa imaginara 12 j   0,sZjeA 1j   11,psPjAe 2j   22psPj Ae AZ   12j X j e  12PA P A De unde: ZA                 ArgXX j j12; 12PAPA Modulul şi faza transformatei Fourier, pentru cazul general, se calculează cu relaţiile: M(1 21) MNkAZ 1k X j k AKrg;             Nkk  11kk25 kAP 1k Transformarea Laplace unilaterala Studiul sistemelor cauzale ce nu sunt initial in starea de repaus (nu au conditii initiale nule) se face utilizand transformarea definita prin relatia :  st          dtetxstx  uL  0 Ea exista si este analitica in semiplanul Re{s}> daca x(t) 0 satisface conditiile : i) Cauzalitate x(t)(t)=x(t); ii) Continue cu exceptia eventual a unei multimi numarabile de puncte in care au discontinuitati de speta intai; iii)  t 0  ,Metx 0,M 0, 0 26 13 Transformarea inversa este definita de: j 1 st1- sRe,dsesXsXtxL 0u 2j j pereche Laplace unilaterala: uL sXtx 27 3 1 Relatia dintre transformata Laplace bilaterala si transformata Laplace unilaterala -Pentru cazul semnalelor cauzale, x(t)=x(t)(t), intre cele doua transformate nu exista diferente :  st st x t s x t e dt x t t e dt  L       st x t e dt x t sL   u  0 -Pentru semnale cu suportul plasat si in t 0, caz în care: 1 , Re 0 ts  L(1 29)   s Dacă treapta unitară este reflectată x(t)=(-t), atunci: 0 1 st st , Re s 0 t e dt e dtt   L    s   şi, deci: 1 L(1 30)   , Re s 0 t        s 1 L L L L L(1 31)    1, , 0t t t                         u u u31 s 5 Proprietati ale celor doua tipuri de transformate Laplace Se noteaza: uLL      sXtx        ssXtx ;DC x LuL     sYty        ssYty ;DC y 1 Liniaritatea:             DC DCx cel putinyax t by t aX s bY s s (1 32)  sbYsaXtbytax uu 2 Translatia in timp: 0st        0 x 0 DC , ,x t t e X s s t st(1 33) intarzierii:  00tsXettxTeorema 00u 3 Modularea in timp: 0st        e x t X s s s00 DC deplasat cu (1 34) ts  0ssXtxe 0u 32 16 4 Scalarea variabilei timp: *1ss , , DC x at X a     a a a (1 35) 1s   0,aXatx  u aa  5 Teorema convolutiei:     DC DCx cel putin yx t y t X s Y s s (1 36)  sYsXtytx uu 6 Derivarare in timp: dx t   DC cel putin sX s s   dt (1 37) ()tdx + ()() 0-↔ xssX u dt33 7 Derivarare in domeniul variabilei s: dX s    s DC tx t    ds (1 38)  sdX  u    ttx ds 8 Integrarea semnalului in domeniul timp: Fie y(t): t  ttxdxty   Se deduce: 1      sXsY s Cum (1/s) are DC Re{s}>0 rezulta: t Xs  s s      DC Re 0 cel putinx d    s  1(1 39)  xsXt0 sXt u u dx saudx   s 00s daca exista un impuls in origine x(-1) (0+) - valorarea integralei in 34 origine 17 9 Teorema valorii initiale a unui semnal cauzal: Pentru semnalul cauzal x(t)=x(t)(t), cele doua transformate sunt identice X(s)X(s) u  (1 40) 0 lim limx sX s sX s   u ss  10 Teorema valorii finale a unui semnal cauzal: Semnalul x(t) fiind cauzal: x(t)=x(t)(t), avem: X(s)X(s) u lim lim limx x t sX s sX s   (1 41)         u    00t s s  35 6 Studiul sistemelor liniare si invariante in timp prin intermediul transformarii Laplace 6 1 Functia de transfer a unui SLIT Pentru un SLIT: y(t)=h(t)x(t) Aplicand transformarea Laplace bilaterala rezulta: (1 42) ;sXsHsY  sthsHL x(t) h(t) y(t)=x(t)*h(t) X(s) H(s) Y(s)=X(s)H(s) Fig 1 7a SLIT caracterizat de Fig 1 7b SLIT caracterizat de raspunsul la impuls h(t) functia de sistem H(s) Raspunsul la impuls caracterizeaza complet comportarea unui SLIT in domeniul timp, iar functia H(s), caracterizeaza complet comportarea sa in domeniul complex Denumiri pentru H(s): - functie (de) sistem sau functie de transfer 36 18 Consideram un sistem SLIT stabil Atunci exista transformata sa Fourier F{h(t)}()=L{h(t)}(j) si axa imaginara este inclusa in domeniul de convergenta al functiei sistem H(s) Exemplu: un sistem anticauzal Conditia de convergenta: Re{s} - 0 F{h(t)}=1/(+j) la impuls h(t) al unui sistem cauzal si functia sistem 0Raspunsul H(s) impreuna cu domeniul de convergenta 38 19 6 3 Sisteme liniare si invariante in timp caracterizate prin ecuatii diferentiale liniare cu coeficienti constanti Fie SLIT caracterizat prin ecuatia diferentiala: kk NM d y t d x t  (1 52)  0 0,kkkkaab 00kkdt dt cu conditii initiale nenule Aplicand transformata Laplace bilaterala, rezulta: NM kk bs Y s s X sa  (1 53) kk 00kk Functia de transfer H(s)=Y(s)/X(s): M k ksb  0kNs Hs  N(1 54)  kDs ksa 0k 39 Radacinile ecuatiei: N k (1 55)    ,0  ksasD k 0 dau polii sistemului, s, in numar de N, incluzand si ordinele de pk multiplicitate: Radacinile ecuatiei: M k(1 56) ,0sbsN  k 0k dau zerourile sistemului, s, in numar de M, incluzand si ordinele de 0k multiplicitate: Se presupune ca sistemul este cauzal (h(t)(t)= h(t)), asadar: H(s)=H(s)=L{h(t)}=L{h(t)} uu Teorema valorii initiale a unui semnal cauzal, aplicata raspunsului la impuls (si el cauzal), conduce la: ssN   lim0 lilim ms H s shHs       u(1 57) sss  Ds  40 20 Daca h(0+) este finit, deci nu exista impulsuri Dirac in origine, atunci se poate deduce ca M+1N Tinand seama de forma expresie(1 54), avem: 12MM     1 1 0MMb s b s b s b s 0h lim   1NN sa s a s a s a     1 1 0NN Daca M+1>N, atunci limita nu poate fi finita Ramane, in consecinta: h(0+) N expresia lui h(t) prezinta impulsuri Dirac in origine Exemplu: 21s    th t e tt       Hs1    u   11ss 41 In mod normal numarul polilor unei functii sistem intrece cu cel putin o unitate numarul zerourilor ei Polii unui sistem stabil sunt plasati doar in semiplanul stang, iar zerourile sale pot fi plasate oriunde Sistemele ce au atat polii cat si zerourile plasate in semiplanul stang se numesc sisteme de faza minima Exemplu Fie: 1 t  sHteth u 1s Raspunsul in frecventa al acestui sistem este: 2 ,1/1jH     arctg  ,1/1H si fie sistemul cu raspunsul in frecventa:  1j   0011tt  12h t e e t ,H         0000    011j 1j 0 42 21 Se poate observa ca raspunsul in frecventa al celui de al doilea sistem are acelasi modul ca si sistemul caracterizat de H() Functia sa faza-frecventa este:  22 arctg arctgarctg          0  00 Functia sistem corespunzatoare lui h(t): 0  0s  sH u 10ss 0 1  sH u 1s (){}(){} ,H s H s   minmin>       0uu Fig 1 12 Functie sistem si 0 domeniul ei de convergenta, DC              00 Prin urmare H(s)=1/(1+s) neavand nici un zero in semiplanul drept, u corespunde unui sistem ce introduce un defazaj minim in comparatie cu alte sisteme (cauzale) ce au aceasi caracteristica H() 43 Contributia polilor unui sistem cauzal la raspunsul la impuls al acestuia Functia de sistem H(s) (1 54) pentru un sistem cauzal, u M k ksb  0kNs Hs  uN(1 54)  kDs kas 0k poate fi descompusa in fractii simple si pusa sub forma: N kA (1 58)  ;sH ,sssHA u sspkuk 1kpksspk daca toti polii sunt simpli Forma lui h(t) este atunci: N ts pk (1 59)  teAth   k  1k Deosebim doua cazuri: i) poli reali s si ii) complecsi s pkpk Cum coeficientii a si b in ecuatia (1 54) sunt reali, rezulta ca polii kk complecsi nu pot aparea decat in perechi complex conjugate 44 22 Poli simpli NN stA k ; pkH s h t A e t            uk    11kkpkss  i) Poli reali, s=  : contributia este o exponentiala pkk kt kAe ii) Poli complex conjugati : s=  j pkkk t sinkB e t    k k k  Poli multipli de ordinul doi A B C D E F          2 2 2* ur i iHs *r i i   p p ps s s s s ss s s s s s         p p p i)Poli reali de ordinul doi: ttr : pp   01pps M e M te ii)Perechea de poli complex conjugati *i i t t                   0 0 1 1, : sin sinpps j s j N e t N te t 45 Pentru un sistem cauzal, stabilitatea stricta este intodeauna asigurata daca polii sunt in semiplanul stang Poli simpli pe axa imaginara = sistem stabil la limita Pentru o singura pereche de poli complex conjugati plasati pe axa imaginara se obtine un oscilator sinusoidal Plasarea unui pol al unui sistem cauzal in semiplanul drept, chiar daca este simplu da nastere unui sistem instabil Polii plasati in semiplanul stang dau o componenta tranzitorie ce se amortizeaza Raspunsul permanent, ce se mentine cu aceleasi caracteristici in timp este contributia exclusiva a polilor simpli situati pe axa imaginara 46 23 Calculul raspunsului unui si sistem caracterizat printr-o ecuatie diferentiala kk txdtydMN aba,0  0kk kk dtdtkk 00 Exemplu: Se considera un sistem de tip cauzal, caracterizat de ecuatia diferentiala: 2 tdytyd ,23txty 2 dtdt Avand conditiile initiale y(0+)=2 ; y’(0+)=-2 Fie x(t)=4(t) Se considera transformata unilaterala:  ;sYty   0'yssYty u u Prin derivare      00yyssYsyy u 47 Transformarea Laplace a semnalului de intrare este: 4   t4  s 2    tdytyd        ,23txty 2 dtdt 4               22322sYssYssYs uuu s Se determina transformarea Laplace unilaterala a semnalului de iesire: 2  ss 222442  sY      2u  ssssss 2123     Se determina raspunsul sistemului: tttt22  teetetetty222222 48 24 tttt22  teetetetty222222 Se pot verifica conditiile initiale: 2tt   2222222limlimeety  00tt Derivata de ordinul intai a raspunsului: 2tt   ,42teety 0  2 tt        eety    24242limlim    tt 00 49 Sisteme de ordinul intai  tdy         txKty(1 66)  0 000 dt K 0        ssHsH(1 67)    Re u 0 s  0 Sistemul:- este stabil, are raspunsul la impuls de tip exponential, polul fiind real Frecventa  este definita de pozitia punctului A pe axa imaginara 1  ;H(1 68) CPZ pentru un PAFig 1 13 sistem de ordinul intai , HPA  Caracteristica de modul este o functie para: PAAP , /2)(/2) -( Caracteristica de faza este o functie impara 50 25 Caracteristici de frecventa K=1 si ω=1 0 ,0→⇒,0=PAA2→ ,1== 0= ⇒∞→/,PA → 0∞ K  ()=0 H j + 0 1 ()= H PA ()↓ ,↑HPA ↑⇒ () -=   ,↑  ↑ (   /2)→ ⇒ ()()/2) -→   ↓ (      1==  PAOA ⇒ ,1== ,2= 4/= 00 51 Hodograf  →   → ,0= ⇒ ∞ ⇒ PA∞ PAA2/-→ → ,∞   ,0→ ,1== 0= 0= 0 1= ⇒1==PAOA ,1== ,2= 4/= 00 52 26 Sisteme de ordinul doi (Studiu individual) 22 tdytydK 22 0 H s H s ,2txKty  000 222u dtdt2ss    00  2 ,1s 0 002,1 2 a) Frecventa de oscilatie:  01 OP 01 Punctul A defineste o frecventa: 1         H; 12 12P A P A         HAA , P si 0P 21   ; si  21 2 >1, cei doi poli sunt plasati pe axa Fig 1 14 a) ξ<1, CPZ Fig 1 14 b) CPZ pentru b)  pentru un sistem de ordinul un sistem de ordinul reala =1, cei doi poli se confunda, in doi avand o pereche de poli doi avand poli reali Daca  complex conjugati punctul - 0 de pe axa reala 53 Cazul polilor complex conjugati Caracteristici de   2/1 si 1 frecventa 0 1        ()()  → ∞ ⇒ → ,∞ → ,∞     ,0=∞ /2,= ,2/=  -=∞     2121HAPAP H; 12 12P A P A ,1==⇒ ()()0=0,1=0 /4,-= ,4/=  0= HAPAP    2121 2-2=⇒ 2+2= ()2=+ tg21=1 tg1,-2=1+2= 1=/,,/,,HAPAP54 212121 27 Hodograf ()()  → ∞ ⇒ → ,∞ → ,∞     ,0=∞ /2,= ,2/=  -=∞ 2121HAPAP0= ,1==⇒ ()()0=0,1=0 /4,-= ,4/=  0= HAPAP    2121  → ∞ 1= ,2-2=⇒ ,2+2= ()2/=+ tg,2/1=1 tg1,-2=,1+2=  1= HAPAP    212121 55 Cazul polilor reali Caracteristici de frecventa Daca =1, apare un pol dublu in - pe axa reala 0  Pentru 10 Tema de curs 56 28 Sistemul “trece tot” js 00  H sHsH u js 00 O sinusoida de orice frecventa, aplicata unui astfel de sistem trece fara ca amplitudinea ei sa fie afectata  trece-tot PZC pentru un sistem trece-tot ZA H     ,1 PA 57            2  OA tgtg OZ 0  arctg 0 Caracteristica de faza () a unui sistem trece-tot Cand =0, deci A0, =0 defazajul introdus de circuit este  Cand OA=,  =/4  ()=-2/4=/2 00 Daca aplicam o sinusoida unui astfel de sistem  nu isi modifica amplitudinea, ci doar faza conform caracteristicii () 58 29 Functia sistem echivalenta unor sisteme conectate in serie si paralel a) Doua sisteme conectate in b) Doua sisteme conectate in serie si sistemul echivalent paralel si sistemul echivalent      ththth   ththth e21 e21             sHsHsH sHsHsH 21e 21e 59 30 